Note on lattice spin in graphene and "spin from 

isospin" phenomenon 

Piotr Kosiriski 1 * Pawel Maslanka 1 } Jadwiga Slawiriska 1 ' 2 * Ilona Zasada 2 
1 Department of Theoretical Physics and Computer Science, 
University of Lodz, Pomorska 149/153, 90-236 Lodz, Poland, 
2 Solid State Physics Department, University of Lodz, 
Pomorska 149/153, 90-236 Lodz, Poland 



Abstract 

It is well-known that the dynamics of low energy electron in graphene honey- 
comb lattice near the K/K' points can be described, in tight-binding approxima- 
tion, by 2+1 massless Dirac equation. Graphene's spin equivalent, "pseudospin" , 
arises from the degeneracy introduced by the honeycomb lattice's two inequiva- 
lent atomic sites per unit cell. Mecklenburg and Regan ( Phys. Rev. Lett. 106 
(2011), 116803) have shown that, contrary to the common view, the pseudospin 
has all attributes of real angular momentum. In some circumstances, the inter- 
nal symmetries can produce an important contribution to angular momentum. 
This phenomenon has been known for many years in particle physics and called 
" spin from isospin" . We show that similar mechanism works in the case of lattice 
pseudospin. 
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Angular momentum is one of the most important notions in both classical and 
quantum physics. This is due to its intimate connection, via Noether theorem, with 
the isotropy of physical space implying rotational invariance of interaction. The angular 
momentum conservation ( together with other invariance principles, like parity) leads 
to the selection rules concerning emission and absorption processes in atomic physics, 
particle scattering etc. Mathematics of rotation group 50(3) ~ SU{2) dictates most 
of the properties of angular momentum: its spectrum, additional rules, uncertainty 
relations for components etc. Another important property of angular momentum is 
that, in general, it is a sum of orbital and spin parts; the former takes only integer 
values while the latter - both integer and half-integer ones. Therefore, once we know 
that the total angular momentum is half-integer it can be taken for granted that the 
spin part gives nonvanishing contribution. 

The separation of angular momentum into orbital and spin parts cannot be inf- 
fered from the structure of 50(3) alone. However, it can be demonstrated by referring 
to the larger symmetry group, Galilean or Poincare one. For example, the structure 
of the unitary representations of the Galilei group (for the Poincare group the situa- 
tion is slightly more involved) implies at once the relevant decomposition of angular 
momentum. 

There is also another way of introducing the spin contribution to total angular mo- 
mentum, called colloquially "spin from isospin" PQ -=- [3] - Some field theories posses 
rich nonperturbative sector separated from the vacuum by a barier of topological ori- 
gin. They are usually characterized by the existence of large internal symmetry. On 
the quasiclassical level nonperturbative sectors are characterized by the existence of 
nontrivial solutions to classical equations of motion. These solutions are stabilized due 
to the existence of the lower bound for energy expressed in terms of some topological in- 
variants. The semiclassical theory is obtained by quantizing small oscillations around 
the classical solution. The symmetry group of the Lagrangian includes 577(2) x G, 
where G is the group of internal symmetries. However, due to the appearance of non- 
trivial soliton solution the actual symmetry group is smaller. For example, if G contains 
577(2) it may appear that the soliton solution is invariant under the 577(2) diagonal 
subgroup generated by the sum of rotation generators and internal SU(2) ones. In 
such a way we obtain additional contribution to angular momentum. As an example 
consider the gauge theory with isospin group as a gauge group. Assume all particles 
entering the theory are bosons so, superficially, there is no spin one-half contribution 
to angular momentum. Consider the monopole solution as "soliton" one. It is "spheri- 
cally symmetric" in the sense that the rotation can be compensated by the appropriate 
isospin rotation. So there is isospin contribution to rotation generator. Now, if some 
matter bosons have isospin | they give spin | contribution to the angular momentum. 
More specifically [2 J , consider the SU (2) gauge theory consisting of gauge field A® and 
two matter fields Q a (isospin one) and U (isospin one half). The model has a magnetic 



2 




Figure 1: Honeycomb crystal lattice of graphene. 
monopole solution of the form 

Q a d = * a Q(r) 

u d = o 

A% = e iab x b W(r) (1) 

the upper indices are the isospin ones while the lower refer to space-time transformation 
properties. So we are dealing with one vector isospin— ~ field, one spinless isospin— 1 
and one spinless isospin— | bosonic fields; no fermions enter the Lagrangian. 

It is now easy to build the quantum theory of small fluctuations around classical 
solutions ([1]). Let A q , Q q and U q be the quantum fluctuations. To lowest order the 
fields A q and Q q do not mix with U q ; in the nonrelativistic limit one classifies the states 
of U q by solving the Schrodinger equation 

HU q = EU q 

H = ^—[Pj ~ eA a clj T a } 2 + binding potential (2) 

It is obvious that A c ij is invariant under simultaneous rotation in ordinary and 
isospin spaces. Therefore, the rotational invariance in generalized sense implies that 
the generator commuting with H reads 

J = xxp + f (3) 

We see that J gets half-integer contribution through T ( T are generators in the 
isospin representation spanned by U so they corresponds to half-integer representation 
of SU(2)). Consider now the electron motion in honeycomb lattice (see Fig.l) 



3 



We are interested in nearest-neighbor approximation. It is easy to write out the 
relevant Hamiltonian [3]. Its main ingredient is the single particle Hamiltonian 



H = -t 
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(4) 



where t is the hopping parameter while A-the energy difference between the sites "A" 
and "5". 

We are interested in the low-energy limit around the K points of Brillouin zone, 
K K = K 2b2 ^ bl + mb\ + n&2, k = ±1. Expanding around these points we find [3]: 
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where k = Q — K. 

Let us note the analogy with "spin from isospin" case. As previously, we are dealing 
with rotationally invariant theory with no half-integer angular momentum (if we neglect 
the processes with spin-flip). The lattice plays here the role of classical background 
configuration which breaks the rotational invariance (due to the appearance of a s and 
a<i in the Hamiltonian (J5])). However, as in the previous case, the invariance can be 
quite easily restored. To see this note the identity 



where 



k^ ■ (ndd — ia s ) = k ■ (na d — ia s )e lK ' p 



k v = (kicoscp + k 2 sin(p, —kisintp + fc 2 cos</?) 



(6) 
(7) 



is the rotated momentum. Now, eqs. (jSJ), © and ([7]) imply 

U(k) = U + H(k v )U((p) 

where 



U{<p) 
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(9) 



Note, that U(ip) plays the role of compensating "isospin" transformation. Eq.® leads 
to the conservation of "total angular momentum" 



J = X X p + K 



f3 



(10) 



Concluding, we see that the lattice pseudospin emerges because the quantization of 
electronic motion is performed on nontrivial background. There is a close analogy with 
the "spin from isospin" phenomena where the half-integer spin results from quantization 
on monopole background. It should be also noted that the above picture is valid within 
the low (electron) energy approximation. If higher order terms of the expansion of the 
Hamiltonian in the powers of ka are taken into account the rotation invariance is broken 
down to the discrete one (reflecting the geometry of the lattice) in a way which prevents 
its restoration by adding some extra terms (like lattice pseudospin). 
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